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Abstract 

Additive manufacturing (AM) could be a novel method of fabricating compos¬ 
ite and porous materials having various effective performances based on mecha¬ 
nisms of their internal geometries. Materials fabricated by AM could rapidly be 
used in industrial application since they could easily be embedded in the target 
part employing the same AM process used for the bulk material. Further¬ 
more, multi-material AM has greater potential than usual single-material AM 
in producing materials with effective properties. Negative thermal expansion is 
a representative effective material property realized by designing a composite 
made of two materials with different coefficients of thermal expansion. In this 
study, we developed a porous composite having planar negative thermal expan¬ 
sion by employing multi-material photopolymer AM. After measurement of the 
physical properties of bulk photopolymers, the internal geometry was designed 
by topology optimization, which is the most effective structural optimization in 
terms of both minimizing thermal stress and maximizing stiffness. The designed 
structure was converted to a three-dimensional STL model, which is a native 
digital format of AM, and assembled as a test piece. The thermal expansions 
of the specimens were measured using a laser scanning dilatometer. The test 
pieces clearly showed negative thermal expansion around room temperature. 
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Additive manufacturing (AM) or rapid prototyping technology has received 
much attention as an innovative manufacturing technique allowing quick fab¬ 
rication of detailed and complicated objects mm- Recent progress of such 
technology has improved the manufacturing accuracy and level of detail to a 
composite or porous material scale (about lO/rm) of a three-dimensional (3D) 
structure. Thus, AM could be applied to the formation of not only macro¬ 
scale parts but also new composite or porous materials. Such materials would 
have various effective properties according to the mechanisms of their internal 
geometries. 

Innovative materials have indeed been developed by integrating mathemati¬ 
cal and numerical design methods and AM technologies; e.g., tissue engineering 
a bone scaffold having compatibility with human bone in terms of stiffness and 
permeability mm and developing materials with a negative Poisson’s ratio mm- 
The greatest advantage of the additive manufacture of innovative materials is 
fast practical realization. In the conventional development of materials, even 
though innovative techniques have been developed in the laboratory, issues re¬ 
lating to bulk production and machining need to be overcome to allow practical 
use. In contrast, a geometry-based material fabricated employing commercial 
AM technology could be directly used in engineering parts. Moreover, function¬ 
ally graduated characteristics could easily be achieved by spatially varying the 
internal geometry [7]. 

A further advantage of AM is multi-material fabrication although this func¬ 
tion is only realized by photopolymer-type AM using commercially available de¬ 
vices. Using multiple materials in the internal structure, there are more degrees 
of freedoms of realizable effective material properties than in the single-material 
case [Sli. Effective thermal expansion is a representative example. By forming 
an internal structure using more than two materials with different positive coef¬ 
ficients of thermal expansion (CTEs) and voids, even negative effective thermal 
expansion could be achieved. With a rise in temperature, usual thermal expan¬ 
sion occurs within the internal structure of materials. However, the mechanism 
of the internal structure converts this thermal expansion into an inward defor¬ 
mation of the outer shell and macroscopic negative thermal expansion is then 
observed. 

The studies of such internal geometry-based negative thermal expansion ma¬ 
terials had been performed by both theoretically [iniiiiiiiiiis] and experimen¬ 
tally [m [THl [H] . By choosing two general materials with known stiffness and 
GTE as base materials, the GTE and stiffness of the composite can be tuned by 
changing the microstructural shape according to elastic mechanics; compounds 
having a negative GTE are limited to special cases and the characteristic tuning 
is an important research topic [niiiH]. However, manufacturability is a major 
issue relating to internal geometry-based negative-GTE materials. First, since 
the internal structure must be composed of at least two materials having certain 
shapes, machining and assembly processes are required for each cell. This might 
be critical in mass production. Moreover, special techniques such as microfab- 
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rication by co-extrusion and reduction sintering m or microelectromechanical 
system fabrication [16] are required to form a small-scale (< 10mm) internal 
structure. 

AM can resolve the above issue since it has the potential to realize multi¬ 
material production in almost the same process as that used for a single ma¬ 
terial by changing the added material according to position. However, only 
photopolymer-type AM is currently available commercially, and the materials 
of such AM have a smaller CTE gap than the combination of metals introduced 
in previous experiments [H nulls]. Thus, techniques are required to design an 
internal geometry that provides effective negative thermal expansion through 
the combination of materials having a small CTE gap. Numerical structural 
optimizations could be powerful tools for such difficult structural design. These 
optimizations could automatically lead to the optimal structure employing nu¬ 
merical structural analysis and optimization techniques. In particular, topology 
optimization could achieve the fundamental optimization of the target structure 
including a change in the number of holes m- Theoretical studies of inter¬ 
nal geometry-based negative thermal expansion based on topology optimization 
have been conducted nniiTs]. 

Against the background described above, the present study develops an in¬ 
ternal geometry-based negative-thermal-expansion material fabricated by mul¬ 
tilateral photopolymer AM. After measuring physical properties of photopoly¬ 
mers, the internal geometry is designed to maximize negative thermal expansion 
employing numerical topology optimization. The effective physical properties 
of the internal structure are calculated employing numerical homogenization. 
Test pieces composed of the designed internal structure are fabricated by pho¬ 
topolymer AM. Negative thermal expansion of the internal structure is then 
experimentally verified by measuring the thermal deformation of the test pieces 
using a laser scanning dilatometer. 

The present work investigates the thermal expansion phenomena of a com¬ 
posite porous material composed of an internal structure that has a periodic 
layout in a plane. We assume the thermal expansion behavior of the material 
in the internal structure follows the constitutive relation 

O'ij — ^ ( 1 ) 

where cr, C, e, a and (3 are the stress tensor, the elastic tensor, the strain tensor, 
the CTE tensor and the thermal stress tensor and AT is a temperature change 
from a reference temperature. In this research, we realize negative thermal 
expansion by considering the effective (or macroscopic, averaging) stiffness and 
CTE of the porous composite. The macroscopic effective physical properties of 
the periodic structure can be calculated employing numerical homogenization 
|20l[2T]. The effective elastic tensor , CTE tensor and thermal stress 
tensor f3^ of the periodic structure composed of a unit cell Y are calculated as 


3 


(2) 



where x Etnd i/j are the displacements obtained by solving the problem oi Y- 
periodic cells expressed as 

= ( 4 , 

(d«-c«„-^)g^<ir = o, (5) 

where v is an arbitrary test function. 

The composite porous internal geometry is designed by optimizing the above 
effective stiffness and CTE. We set a square domain as the base shape of 
the internal geometry and design it by allocating two types of photopolymers 
and voids. Solid isotropic material with penalization (SIMP)-based multiphase 
topology optimization is introduced as the design optimization tool for the inter¬ 
nal geometry, which is represented by the layout of three phases, two materials 
and voids, in the specified domain by defining two artificial density functions </>! 
and (j )2 (0 < (^ 1 , <^2 < !)• The elastic modulus and CTE are formulated as func¬ 
tions of (^1 and 4 > 2 i which represent the existence of materials and the kinds of 
materials respectively. The local Young modulus and CTE of the design target 
domain are represented as 


E{(t)l,4>2) = 4>\{4>2 Ei + {I - (j)2)E2} , ( 6 ) 

q;((/> 2) = (^2ai-I-(1 - </>2)a2, (7) 

where Ei, E 2 , oi and 02 are Young’s modulus and the CTE of photopolymers 1 
and 2 respectively. In other words, (pi k. \ and (/)2 ~ 1 means that photopolymer 
1 exists, (pi K, 1 and p 2 ~ ^ means that photopolymer 2 exists, and ~ 0 means 
the void exists. Since the state cannot be identified from intermediate values 
between 0 and 1, it should be avoided by contriving the problem and optimizer 
settings. To realize the isotropic effective CTE, symmetry is assumed on the 
center and diagonal lines of the square design domain. Thus, distributions of 
pi and p 2 are optimized only on the 1/8 domain shown in Fig. [5] (a). 

When only the effective CTE is considered in designing a structure having 
negative thermal expansion, a structure having very low stiffness could result, 
which would be unsuitable for fabrication and experiment. Moreover, the effec¬ 
tive thermal stress tensor (3^ is preferred over the effective CTE both in ([3]) 
in terms of the easiness of the calculation of its gradient used in the optimizer. 
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We thus intend to maintain a certain level of effective stiffness while reducing 
the effective CTE by setting the multi-objective function as 

minirnize = -w * + {1 - w) * (i = 1, 2). (8) 

9l,<p2 

The steps to the optimization are as follows. We first solve O and 
employing the finite element method (FEM) and a commercial FEM solver 
COMSOL Multiphysics (COMSOL Inc., USA). Second, the effective physical 
properties in ([2]) and ® and the objective function in ([8]) are calculated. Since 
the density function is updated by gradient-based algorithms, the first-order 
gradient of the objective function is then calculated using the adjoint method 
[22] . The density functions (f>i and ^2 are updated by the method of moving 
asymptotes (MMA) [531 stage of the optimization. In the latter 

stage, to obtain a clear shape avoiding intermediate values, the density function 
is updated using the phase field method [24] . 

According to the numerical design methodology described above, the pho¬ 
topolymer composite internal geometry is designed. An Objet Connex 500 
(Stratasys Ltd., USA), which is the only commercial photopolymerization man¬ 
ufacturing machine offering multi-material 3D printing, is used to fabricate the 
test pieces. The machine produces 3D structures by spraying liquid photopoly¬ 
mer onto a build tray in thin layers and exposing the photopolymer to ultravi¬ 
olet light. In producing a structure from multiple materials, a rigid material, a 
rubber-like material and several admixture materials can be used. 

We measured Young’s modulus by tensile testing in a temperature-controlled 
bath at 20, 30 and 40°C and measured CTE using a connecting rod dilatometer 
for several materials between RT and 50°C, which is close to 45°C, the heat 
deflection temperature (ASTM D648, 0.45MPa) of a rigid material. Accord¬ 
ing to the obtained physical properties, we chose VeroWhitePlus RGD835 and 
FLX9895-DM, which is an admixture material of VeroWhitePlus and TangoB- 
lackPlus, because this combination achieves a certain level of stiffness and a 
certain CTE difference in the temperature range between RT and 40° C. The 
measured physical properties are listed in Tableland plotted in Fig. [T] Pois¬ 
son’s ratios of these materials were measured only at 40° C and they are both 
about 0.48. Although CTEs of both materials have strong temperature de¬ 
pendency and dispersion within this temperature range, we found a difference 
in the CTE on average, which is a necessary condition for a negative effective 
CTE. With the aim of realizing a negative effective CTE at 40°C, which is the 
temperature at which the CTE difference between the materials is a maximum, 
we set Young’s modulus and the CTE in ([6|) and o as follows according to the 
measured average physical properties: 

El = 290MPa, E 2 = 5.0MPa,ai = l.OxlO-'^K"^ and 02 = I.2x10“'‘K“^ (9) 

Topology optimization is then performed using the optimal shape shown in 
Fig. 2(a) of ref[T0] as the initial shape. We first set the physical properties as, 
El = 290MPa, E 2 = 290MPa, oi = 1 x and 02 = 1 x which 
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Table 1: Measured Young’s modulus of VeroWhitePlus RGD835 and FLX9895- 
DM photopolymers. iV = 3 for each temperature. 


Temperature [°C] 

VeroWhitePlus RGD835 
Average [MPa] SD [MPa] 

FLX9895-DM 
Average [MPa] 

SD [MPa] 

20 

1996.2 

111.2 

67.8 

2.4 

30 

1023.3 

228.0 

13.0 

5.4 

40 

288.5 

34.6 

3.8 

0.8 


Table 2: Effective planar material properties of the optimal internal structure 
with the bulk physical properties in (jH]) 


Effective elastic tensor [MPa] 


Effective GTE tensor [K ^] 

C’iiii(= C2222) C'f^22(= ^2211) 

^1212 

«fl(= «2^2) 

4.12 X 10“^ -1.01 X 10“^ 

1.90 X 10“^ 

-3.53 X 10-'’ 


has the same ratio of each E and a with the example used in ref [10]. They 
were varied to the values in ([9]) gradually during the optimization procedure. 
Figure [2](a) shows the final converged solution. We had introduced the forced 
void domain (marked by gray dots) in which no material can exist in the final 
stage of optimization to maintain enough space for thermal deformation. 

We then obtained 3D data in STL format by extruding the two-dimensional 
(2D) optimal shape in the thickness direction. The 2D view of the STL is shown 
in Fig. [2](b). The hinge-like thin rigid photopolymer 1 on the edge of the diag¬ 
onal seemed not to endure the fabrication and was replaced by the thicker soft 
photopolymer 2. The effective in-plane elastic tensor and GTE tensor calculated 
for the 3D model with the physical properties in (|9|) are presented in Table [2] 
Negative effective GTEs were certainly obtained. An effective negative Poisson’s 
ratio was also incidentally observed. The free thermal deformation shape of this 
model is shown in Fig. [2](c). The mechanism of the negative thermal expansion 
is as follows. Bending deformations occur in the bilayer composed of different 
GTE materials. Since these parts are connected by a hinge-like part, the parts 
connecting with the next cell deform inward. When the inward deformation of 
the connecting parts is larger than the axial thermal expansion of the connect¬ 
ing parts, there is macroscopic negative thermal expansion. According to the 
homogenization analysis for the STL model using Young’s modulus shown in 
Table [1] and the parameterized GTE, differences in the GTEs of photopolymers 
of at least 1.15 times is required at 20, 30 and 40 °G to realize the negative 
effective GTE. 

We employed laser scanning to measure the thermal deformation since it 
is suitable for relatively large specimens. The testing device was an SL-1600A 
(Shinagawa Refractories Co., Ltd.), which measures the axial thermal deforma¬ 
tion of rod-shaped specimens by laser scanning. 

By arranging the eight base STLs shown in Fig. [2] (b) in two lines, the STL 
model of the test piece was constructed. The size of each internal structure was 
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10 X lOtzmeslOmm and the total size of the test piece was 40 x 20timesl0min. 
Figure [3] shows the outline and close-up pictures of the fabricated test piece. 
From the close-up picture of the internal geometry, it can be said that almost 
the same shape as the original STL model was achieved except for small branch¬ 
like features. 

The thermal expansions of the long sides of three test pieces were measured 
in the temperature range between RT and 50° C with a temperature rising rate 
of l°C/min. The measured strains are plotted in Fig. 2] Negative thermal 
expansion was clearly observed for each result between RT and about 45°C. The 
early thermal expansions are linear, and corresponding CTEs were calculated as 
1.12-1.18K“^ between RT and 34°C. Since the variation in CTEs is relatively 
low in such a temperature range according to the bulk CTE results shown in 
Fig. m stable differences of CTEs of the two materials are guessed for this 
linear region. However, beyond approximately 34 to 36°, the CTE differences 
decreased and negative thermal expansion reduced or became positive. Owing 
to the variability of the physical properties of photopolymers, the performances 
of the fabricated porous composite also varied widely. However, planar negative 
thermal expansion characteristics could be clearly designed for photopolymer 
porous composites. 

In summary, we fabricated a porous material with planar negative thermal 
expansion employing multi-material photopolymer AM. The internal geometry 
was designed using topology optimization. The mechanism of the effective neg¬ 
ative thermal expansion is inward deformation resulting from the combination 
of bending and hinge mechanisms. The 2D optimal internal geometry was con¬ 
verted to an STL model and assembled as a test piece. The thermal expansion 
of the test piece was measured with a laser scanning dilatometer. Negative 
thermal expansion corresponding to less than —1 x was certainly ob¬ 

served for each test piece of the fV = 3 experiment. Owing to the temperature 
dependence and variable physical properties of the photopolymer, the control 
of the performance of the developed materials appears difficult. However, al¬ 
though multi-material AM is currently commercially available only in the case 
of a photopolymer, multi-material metal AM is under development [151 [55]. If 
such AM is realized, the design process proposed in this paper could be used 
for materials that are more stiff and stable than photopolymer. 

This work was supported by JSPS KAKENHI Grant number 15K12557. 
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Figure 1: Measured CTEs of bulk VeroWhitePlus RGD835 and FLX9895-DM 
photopolymers. The temperature rising rate is l°C/min. The temperature 
range is between RT and 50°C. N = 3. 
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Figure 2: Optimal configurations and 2D view of the STL model, (a) Original 
topology optimization result. Pink dotted lines indicate the optimized domain. 
Gray dotted lines indicate the forced void domain for maintaining the defor¬ 
mation range, (b) 2D view of the STL model generated from the topology 
optimization result. The hinge-like part of the stiff material on the corner was 
replaced with thicker soft material as shown in blue dotted circles, (c) Thermal 
deformation shape obtained from the reanalysis of the STL model. The green 
lines indicate the original shape. 
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Figure 3: Picture of a test piece, (a) Outline view, (b) Close-up view corre¬ 
sponding to the STL model shown in Fig. [2] (b). 


y = 2.04x10'^ - 1.13xl0“\ R^= 9.90x10 * 
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Figure 4: Thermal deformation of test pieces measured by a laser scanning 
dilatometer. The temperature rising rate is l°C/min. The temperature range 
is between RT and 50°C. RTs of measurements No. 1, 2 and 3 are 17.45, 19.45 
and 22.85°C respectively. Approximation lines are plotted for the data within 
the temperature range between RT and 34° C. 
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